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ABSTRACT 

Meridional flow is thought to play a very important role in the dynamics of the solar 
convection zone; however, because of its relatively small amplitude, precisely measuring 
it poses a significant challenge. Here we present a complete time-distance helioseismic 
analysis of about two years of ground-based GONG Doppler data to retrieve the merid¬ 
ional circulation profile for modest latitudes, in an attempt to corroborate results from 
other studies. We use an empirical correction to the travel times due to an unknown 
center-to-limb systematic effect. The helioseismic inversion procedure is first tested and 
reasonably validated on artificial data from a large-scale numerical simulation, followed 
by a test to broadly recover the solar differential rotation found from global seismol¬ 
ogy. From GONG data, we measure poleward photospheric flows at all latitudes with 
properties that are comparable with earlier studies, and a shallow equatorward flow 
about 65 Mm beneath the surface, in agreement with recent findings from HMI data. 
No strong evidence of multiple circulation cells in depth nor latitude is found, yet the 
whole phase space has not yet been explored. Tests of mass flux conservation are then 
carried out on the inferred GONG and HMI flows and compared to a fiducial numerical 
baseline from models, and we find that the continuity equation is poorly satisfied. While 
the two disparate data sets do give similar results for about the outer 15% of the interior 
radius, the total inverted circulation pattern appears to be unphysical in terms of mass 
conservation when interpreted over modest time scales. We can likely attribute this to 
both the influence of realization noise and subtle effects in the data and measurement 
procedure. 

Subject headings: Sun: atmosphere — Sun: interior — Sun: oscillations 
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Introduction 


It is well established that an accurate measurement of solar meridional circulation would be 
of enormous value as an input t o the complex mod els used to study the solar dynamo and large- 
scale convection-zone dynamics (jCharbonneaull2014l l . Furthermore, debate still continues as to the 
precise mechanism for generating and sustaining these poleward a nd equatorward flows and their in¬ 
tricate interplay wit h differential rotation and thermal structure (lRempelll2nnfil : iMiesch fc Hindman 


2011 


DikDatill2014l ): different parametrizations in numerical models produce drastically different 


outcomes. A better understa nding would also have significant impacts for studie s of other stars 
with surface convection zones (jBrun fc Palacios! 120091 : iFeatherstone fc Mieschl 120151 ). 


Over the past two decades, many observational attempts have been made to quantify this 
important phenomenon. Methods us ed have generally included Doppler-shift diagnostics and 


2012 

) for relatively near-surface studies, the various tools associate 

d with local helioseismology 

dGiles et al.lll997: 

Braun & Fanil 1998: 

Gonzalez Hernandez et al.lll999 

Woodard et al.ll2012: 

Komm et al. 

20ld: 

Zhao et al. 

201.31. 

2014 

), and even global helioseismology ( 

Schad et al. 

2011. 201,31 to probe 


the entirety of the convection zone. While a detailed comparison of all of the results from these 
analyses will not be presented here, it is reasonable to state that no consistent, coherent, nor satisyf- 
ing picture of the full circulation prohle has yet to emerge. In particular, observational signatures 
that distinguish between the canonical single, large meridional cells or multiple cells stacked in 
radius or latitude have not become readily apparent. 


In a recent paper (jKholikov et al.l 120141 . hereafter Paper I), measurements of travel-time dif¬ 
ferences of almost two years of ground-based Global Oscillation Network Group (GONG) data 
showed evidence of the signature of a rather complex structure of the large-scale flows in the solar 
convection zone. Here, we add to this rich body of literature and introduce a helioseismic inversion 
procedure, test it with numerical simulations and “known” results, and hnally invert the GONG 
travel times to recover flows. Without any strong reason to think that our results should signif¬ 
icantly differ from others, or will convincingly demonstrate once and for all what the real Sun 
is doing, we nonethe less find reasonable agreement with th e space-base d Helioseismic an d Mag¬ 
netic Imager (HMI, IScherrer et al.l 120121 : ISchou et al.l 120121 ) inferences (jZhao et al.l l2010l ) in the 
near-surface after applying a similar systematic correction to the travel times. 

§ [2] briefly summarizes the GONG measurements and § [3] describes the forward and inverse 
method we apply. We test the inversion technique in § 0] with artificial data from a numerical 
model of meridional circulation and using the well-known differential rotation profile. While we 
do not discuss at any great length how the results shown here fit into context of the overall body 
of literature on this subject mentioned above (this is reserved for an upcoming publication), in 
§ El we present our results, but then seriously question how realistic they are with regards to mass 
conservation in the convection zone. Gonclusions are provided in § [6l 
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Data and Measurements 


We use the measurements that were described in Paper I. To briefly summarize, 652 days of 
tracked GONG Doppler velocity images wer e used to construct t ravel-time differences using the 
spherical harmonic decomposition technique (jKholikov et alJl201ll ). Skip distances from 3° to 47° 
were considered, which correspond approximately to lower turning points of 0.98 to 0.70 
(angular degree 0 — 300). After averaging over longitude, the final “coverage” of the south minus 
north (S-N) travel times on the solar disk is within the latitude range of about ±75° for the shortest 
distances to about ±40° latitude for the longest. 


To correct for the center-to-limb (CTL) systematic effect as first presented in IZhao et al 


(j2012al i. the identical measurement procedure was carried out to obtain east minus west (E-W) 
travel-time differences. Essentially, we simply rotate the Doppler images by 90 degrees and use the 
same data-analysis pipeline. 

Results were shown in Figures 2 and 3 of Paper I. What we found was that, unsurprisingly, 
initial S-N travel-time differences gave counterintuitive results, in that the signal became stronger 
at larger depths (skip distances). Presumably, as mentioned above, this is due to an unknown 
artifact that trends in the center-to-limb direction. After the subtraction of the E-W measurements 


from the S-N ones, which has been argued will correct for such a trend (jZhao et al.ll2012al i. the S-N 
travel times show interesting behavior that suggests a signature of the equatorward return flow, and 
possibly multiple cell structures. In the inversion analysis presented here we use the E-W corrected 
S-N travel-time differences from Paper I, and revisit the systematic effect in the discussion section. 


3. Forward and Inverse Methods 


3.1. Sensitivity Kernels 


We work in the acoustic ray theory (e.g. iKosovichev &: Duvalll Il997l i. whereby travel-time 
differences 5t of waves propagating in opposite directions are assumed to be caused only by slowly- 
varying interior flows along a given ray path. The formalism is usually expressed as 


5t = —2 


f n - u 

Jr c| 


ds, 


( 1 ) 


where u is the flow, n is a unit vector tangent to the unperturbed ray path T and Cg is a model 
solar sound speed profile. Working in a coordinate system where r is the distance from the Sun’s 
center and 6 the latitude, Eq. ([T|) can be recast in the form 


6t = J K{9,r) ■ u{9,r) ds, (2) 

where K are known as sensitivity kernels. For the study here, only the horizontal, or meridional, 
component of the kernels and flows is considered. Ignoring the radial component of the flow is 
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presumably justified since its influence on the travel-time differences is small, although further 
discussion of this is given in §[5j To interpret the GONG travel times, we compute kernels for each 
measurement, namely for all 45 values of A in the range of 3°-47°, using 3300 ^uHz for the value for 
the central frequency. 


3.2. SOLA Inversions 

We solve Eq. ([2]) by inverting the travel-time differences and obtain the meridional flow as 
a function of latitude and depth using the subtr active optimally localised averaging, or SOLA 
(jPiipers &: Thompson!Il994l : iJackiewicz et al.ll2012l l technique. To simplify the description we in¬ 
troduce the following notation: let be the set of measured travel-time differences, a = {Aobs} 
the set of skip distances, and j3 = {0obs} the set of central latitudes. 

Let the position of interest at which we would like to estimate the flow speed u in the solar 
interior, i.e., the target location at some depth and latitude, be The goal of the inversion 

is to And a linear combination of the travel times 

u{9T,rT) = Wa^i3{9T, rT)ST°‘’^ 

a^l3 

with yet to be determined weight functions Wa,j 3 - To derive the weights, consider the cost function 




K°‘'^{9, r)wa,/3{9T, rT) — T{9, r; 9t, rx)) 


a,/3 


d9 dr 


Wa'^l3^{0T^rT), (3) 

a,/3 a',13' 

where are the sensitivity kernels described earlier and // is a trade-off parameter that con¬ 

trols localization and error propagation. T(0,r; 0 Tj^t) is the target function, typically a Gaussian 
centered at the position of interest that is normalized as 

J j T(0, r; 0 xTt) d0dr = 1. 

Averaging kernels, seen as the first term in Eq. ([3|), are similarly normalized: 

f f K{9,r;9T,rT)d9dr = f / ^ r) d0 dr = 1. 

Je Jr Je Jr ^ 

a^p 

With these constraints and an added Lagrange multiplier A, the minimization is performed on the 
the functional 


9(u;; A,/r) = x(u;;/r) -f A 


K{9, r; 9t, rT')d9 dr — 1 ) , 


9 Jr 


(4) 
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whereby ones solves the coupled equations 


— - n 

dwaA^T,r^ ’ dx 


(5) 


The travel-time noise covariance matrix, which is used to compute the variance of the noise prop¬ 
agated in the inversion (second term in Eq. [3]) , has the form 


Kfi 3 ' = Cov 




The error in each travel-time measurement, due to random realization noise, is captured in N. In 
our particular analysis here, this covariance matrix is diagonal, and since the kernels and 

measurements dr are normalized by the error in the set of observations, the covariance matrix is 
simply the identity matrix. 


Given this inverse problem the solution can be found as 

Wa,l3{0T,rT) = 

where 

Aa,p = jJ dMr + 

ha,f} = jj^ r)T{e, r; 0t, tt) d9dr. 

When solving the inverse problem with the SOLA method, is often poorly conditioned and 

its inverse can be found using a single value decomposition (SVD) method, for example. There 
are two ways to smooth the resulting solution: (i) use two parameters for regularization: // for 
the error covariance matrix, and some threshold paramenter for the eigenvalues to compute the 
inverted matrix using the SVD algorithm; (ii) increase the values of the target Gaussian sizes in 
the two relevant dimensions, FWHM^ and/or FWHMj.. In this work we used different values for all 
of these parameters to find an optimal set, which is chosen by inspection of the so-called L-curves. 


4. Validation Tests 

We provide two examples to test and validate the procedure described above. In this section 
we apply helioseismology to artificial data from a numerical model of enhanced meridional flow, 
followed by an inference of the well-studied solar differential rotation profile. We hnd that within 
the errors the method acceptably recovers the known flow profiles. 


4.1. Numerical Simulation 


The artificial data is described in lHartleo et al.l (j2013l L This simulation has an imposed merid¬ 
ional profile as shown in the leftmost panel of Fig. [2l It has the “classic” structure of poleward 
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Fig. 1.— Forward modeled travel-t ime differences in th e S-N sense (left) and measured ones (right) 
for the simulated meridional flow (jHartlep et ahlLOla ) as functions of latitude and skip distance. 
The color scale is the same in each panel. 


motions in the outer 15% in radius and an equatorward flow starting at depths of 150 Mm in each 
hemisphere. To increase the signal-to-noise without running the simulation for an unreasonable 
amount of time, the authors increased the model flow speed (from its original formulation) by a 
factor of 3 6 to give a surface m aximum value of 500 ms“^, and a maximum return speed of roughly 
lOOms”^ ( Hartlep et al. 2013 1. 


Travel-time differences in the south-north sense are measured for 198 latitudes in the range 
of —85.695° to -1-85.695° and 29 distances in the range of 4.5° — 46.5°. Flow sensitivity kernels 
are computed for each distance, and are used to compute forward travel times using Eq. ([2]). A 
comparison of the two sets of travel times are given in Fig. [TJ As expected, forward travel times 
are smooth (the noise is not modeled). The measured travel times show various structure due to 
noise, and are furthermore not symmetric about the equator, even though the flow model is. Errors 
are not shown but rms values are on the order of about 1 s at all distances, and are not a strong 
function of latitude since this is not a real so lar data set. These results are comparable to the 
helioseismic measurements that were shown in iHartlep et al.l (j2013l l. 


Both the forward and measured travel times are inverted using the methods described in § 13.21 
to recover the strong flows in the model. An optimal set of weights is computed and used to combine 
each set of travel times. The results of the inversions are shown in Fig. [2l Excellent agreement 
is seen for the case using the forward travel times (middle panel), thus providing a convincing 
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Fig. 2.— Simulation and inversion results. Shown are meridional cuts of the flow model (left), 
the inversion of forward travel-time differences (middle), and inversion of measured travel-time 
differences (right). For reference, dotted radial lines are plotted every 15° and dotted concentric 
lines are shown for depths r = [0.7,0.85, l.O]i? 0 . Not shown is the model radial flow profile. All 
plots use the same color scale, but note the different depth and latitude ranges. 

validation of the inversion method and codes. Due to somewhat noisy travel times, the flows shown 
in the rightmost panel of Fig. [2] exhibit features that are not completely consistent with the model. 
The return flow is detected; however, it is truncated in its latitudinal extent, most noticeably near 
the deep equatorial region. The inversion errors in the region deeper than O.85i?0 have an rms of 
45ms“^, while those shallower than that are about 20ms“^. 

To some degree, and without an obvious alternative on hand, this exercise justifies the ray 
approximation for our purposes in modeling travel-time differences related to large-scale meridional 
circulation. It also allows us to find optimal inversion parameters that will be applied to GONG 
data in §[5l 




















4.2. Recovering Differential Rotation with Time-Distance Helioseismology 


The solar differential rotation obtained from global- mode seismology (frequency splittings) is a 
very robust measurement (jSchou et alJll998l : lHowell2009l l , and serves a s a useful test for other meth¬ 


ods. Local helioseismology has also shown success in this realm fe.g. iGiles et alJll998l : iBasu et al 


19991 ) ■ We attempt a validation of the specific methods described here using real solar data from 


GONG. 


In addition to the E-W measurements discussed in §[2] of tracked data, we also made travel¬ 
time measurements of untracked data that preserve rotation information. The geometry is such 
that cross correlations for a given latitude and separation distance are computed from point-to- 
arcs, where the point is always on the central meridian and the arcs lie in one of the hemispheres 
centered on that latitude. The westward signal (prograde) is subtracted from the eastward signal 
(retrograde), leaving behind the expected influence of rotation that is subsequently free from the 
GTL effect described earlier. In other words, since a purely symmetric GTL effect has the same sign 
in each hemisphere with this measurement geometry, upon subtraction it is canceled out. These 
calculations are then done for many distances and latitudes, for about 360 days of data. 

To invert for solar rotation we use the following set of parameters: the location of the target 
function in depth is set from O.SRq to I.ORq with steps of O.OIRq and in latitude from —60° to 
-|-60° with steps of 2.5°; the width of the target function in depth and latitude is FWHM^ = 0.04Rq 
and FWHMg) = 5°, respectively. To make a more meaningful comparison, the global helioseismic 
results are smoothed to our inversion resolution by convolving with the target function used at each 
depth. Finally, both sets of results are averaged over three latitude bins from 0° — 20°, 20° — 40°, 
and 40° - 60°. 


The com parison with the t ime-averaged rotation profile from global analysis of about 1000 days 
of HMI data (IHowe et al.ll201lh is shown in Fig. O Within the uncertainties of our local-helioseismic 
inversion, the results are in agreement over most depths and latitudes except for the near-surface 
region around the equator. There are many possible contributions to any disagreement in general, 
including different instruments and non-overlapping time series. Despite some discrepancies, this 
provides more confidence that large-scale flows, albeit strong ones in this case, are somewhat 
accurately attainable with these tools. 


5. Results and Discussion 


5.1. Meridional Flows from GONG Data 


The GTL-corrected S-N travel times measured from GONG data and intro duced in ^[2, were 
inverted for flows and are shown in Fig.[31 Also provided are the HMI results from IZhao et al.l (j2013l l 
that were kindly provided by J. Zhao for comparison. There are several noteworthy remarks: 
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Depth [Rq] 


Fig. 3.— Comparison of the solar (differential) rotation velocity estimated from inverting travel¬ 
time differences (dashed lines with circles) with global helioseismic results (solid lines). Averages are 
computed over three latitude bins as labeled. Inversion errors in inferred velocity and uncertainties 
from the target function width are shown every three depth points for clarity. 


• GONG and HMI show very good agreement regarding the location of the change of sign from 
surface poleward flows to equatorward flows at about r = O.QI/Zq. This is consistent in both 
hemispheres, as can be seen in panels A, B, and C of Fig. [H 

• The peak near-surface poleward flows occur at about 30° latitude in each hemisphere for both 
data sets (panel D), with amplitudes in the 15 — 20ms“^ range. 

• The shallow return flow peaks at similar latitudes with amplitude around 5 — 10ms“^. 


Unlike HMI, GONG only shows weak evidence of a second “cell” at low latitudes within 
15° of the equator, a feature we suspect is spurious. The equatorward circulation at the 
bottom of the dete ction region is consistent with recent predictions from flux-transport models 
( Hazra et al. 2014 1. 

Any features in high-latitude regions (> 55°) are suspicious, particular the strong equatorward 
flows in the GONG results. 


The results of HMI and GONG are basically anticorrelated at the deepest layers, yet both 
interestingly cross zero at the same latitude (20°, panel E). 
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Fig. 4.— Comparison of meridional flows measured from GONG and HMI data. Panel A shows 
cross-sectional views of the meridional flow profiles within latitudes of ±70° and depths above 
O.74i?0. The dotted lines on each image are plotted at depths r = (0.76,0.85,0.92, l.O)/?© and at 
latitudes 9 = ±(10°, 40°, 60°), for reference in the line plots on the right. Panels B and C show 
velocities as a function of depth averaged over 10° latitude bands, as labeled in the legend, for 
the nothern and southern hemispheres. Panels D and E show velocities as a function of latitude 
averaged over the depth ranges shown in each legend. Amplitude error and radial resolution bars, 
shown only for GONG every several points, are similar for HMI. 


We do not observe any evidence of multiple cells in latitude at any depths, unlike the 


Schad et al.l (2013j) study. 


There are also several (of likely many) caveats to keep in mind: 


Due to the use of long time series, the random errors propagated through the inversion 
are small. However, we c aution that uncertainties due to known and unknown systematics 
(iDuvall fc Hanasogd 120091 ) are likely much larger, and would particularly affect the results 
deep in the convection zone and at high latitudes. 


• That the near-surface meridional flow is not strongest at the very topmost boundary is due 
to the fact that the smallest skip distance is still rather large (using medium-^ data), and 
that there is a level of smoothing in depth due to averaging kernels in the inversion. 
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• To do a more detailed comparison, several independent sets of measurements are needed, for 
example, analysis of separate four-year data sets from 1996-2014 using GONG, MDI, and 
HMI would be ideal . Efforts are currently underway to do this. 


Given these methods and results, one can certainly argue that the ad-hoc correction of the 
CTL systematic is problematic, as the signal one is searching for in the large-distance (deep) travel¬ 
time differences is of order one second. Subtracting sets of large numbers to reveal t his signal is 


not a r ecommended practice, and even though it was convincingly demonstrated in IZhao et al 


()2ni2bl ) that this correction brought multiple HMI observables into agreement, further justification 
is certainly needed. One could imagine that the efficacy of the systematic effect removal should be 
shown first by testing it to recover well-known features, such as the solar rotation. Unfortunately, 
the geometry of the E-W measurements removes the effect in each hemisphere as mentioned earlier. 


5.2. Mass Conservation in the Convection Zone 


It is important to consider if the resulting steady flow fields from the inversions are physically 
plausible. What is not typically shown in the relevant helioseismic literature is how derived flows 
satisfy our notions of a convection zone where material does not enter or leave to any large degree. 
Although it is p ossible to u se mass conservation as an additional constraint directly in the inversion 
equations (e.g. iGilea 120001 ) . we have not implemented that. At zeroth-order, we have already 
ignored inverting for the radial flow. Here we discuss if it even matters to do detailed comparisons 
of meridional circulation solutions since there may be larger issues to first identify. 


We start with the continuity equation 

dp 

di 


+ V • (pv) = 0 


( 6 ) 


for some velocity held v. Ignoring any latitudinal or azimuthal dependence of the density p, the 
expression for mass conservation in spherical coordinates can be expressed as 


dp 

Vr ■ {pVr) + Ve • {pVg) = - — , 


( 7 ) 


or. 


1 

-Vr 

P 


dp 1 5 . 2 N 
^ + u) + 


dr 


dr 


1 d 

r sin 9 dO 


(sin^ue) 


1 dp 
p dt ’ 


( 8 ) 


where 9 is latitude. Since it is not expected that the Sun’s density prohle is changing over moderate 
time scales, the degree to which the total divergence of the mass hux pv is zero is a good indication 
of satisfying continuity. In what follows we consider applying this equation to hve meridional 
flow prohles. Since three of the prohles described below are from inversions that do not have 
a radial velocity component, we compute it from solving the continuity equation requiring mass 
conservation (to some numerical precision). The solution is found using a simple fourth-order 
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Runge-Kutta technique that initializes the radial flow to be zero at the surface boundary point 
(r ~ I.ORq) and integrates from that point inwards toward the deepest part of the convection zone 
considered. Once determined, we may then ask if this flow is indeed reasonable. 


5.2.1. Fiducial model: MFl 


To test such computations and to establish a numerical baseline, we first utilize a model of 
meridional circulation that by construction analytically satisfies Eg. (lUD. Such models are com¬ 
mon i n the literature, and we choose the one originally described in Ivan Ballegooiien &: Choudhuri 
(j 19881 1. which has the “canonical” single-cell deep return flow, and uses a power law density profile. 
We compute this model, denoted “MFl,” and run it through the continuity equation to determine 
the numerical precision to which it is satisfied. The latitudinal and radial flows are shown in columns 
1 and 2 for model MFl in Fig. [5l The sum of the horizontal and radial mass flux divergences scaled 
by the inverse density, i.e., an estimate of p~^{dp/dt), is shown in the 4th column. 


As a test of the initial-value solver, we also compute the radial flow obtained only from (the 
same) vq. This will be denoted Vr in each case hereafter. Column 3 of Fig. [5] shows this derived 
quantity, and we see that it matches very well the original Vr (column 2). For completeness, we 
then recompute the continuity equation given the new W and original vq, shown in column 5. It is 
easy to see by eye that this model satisfies continuity, as it must by construction. 


This template is then followed for 4 more cases below, where all profiles are computed on (or 
interpolated onto) a grid of 350 points in depth and latitude in the convection zone (experiments 
with more points were done but do not qualitatively change the conclusions). In addition to Fig. [5l 
Table [T] shows the quantitative values for a more detailed view of these calculations. Keeping with 
the example for MFl, the 6th row of the table, which is the sum over the mass-flux divergence, 
gives the numerical precision we anticipate in this best case, ~ 10“® in arbitrary units. The 8th 
row gives the same quantity using the recomputed radial velocity, which in this case is very similar. 
Rows 9-10 are important for this discussion since they show the relative magnitude of the numerical 
integral of the two terms in the equation with respect to one of the terms. 

Row 11 of Tabled] is the inverse of the RHS of the density-scaled continuity equation, which 
gives units of time. One plausible interpretation of this is that it represents the timescale over which 
it would take to drain (positive value, dp/dt < 0) or fill (negative value, dp/dt > 0) a “typical” 
region of the convection zone of that model, with mass. This number is computed from the median 
of all the values (column 4 of Fig. [5|) throughout most of the convective domain. So it represents 
an approximate measure of the lifetime for the bulk region of the convection zone. The MFl model 
sets the baseline for the numerics at about 10 kyr. 


Another appealing way to think about this is if we take the MFl model as ground truth (in 
a numerical sense, by construction) to represent perfect mass conservation in the convection zone 
over about Iq = 4.57Gyr, we find a scaling factor (~ 10^). This scaling can be applied to the 
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subsequent MF profiles to understand the calculations in terms of the solar lifetime. This timescale 
is shown in row 12 for all cases. 


5.2.2. Profiles M2-M5 


We next con sider thelHartleo et al.l (1201311 flow model, “MF2,” employing the standard Model 
S density profile (jChristensen-Dalsgaard et al.l 119961 ) for the computations. As can be seen in the 
4th column of Fig. 0 for MF2, the mass is conserved very poorly. We speculate that this is somehow 
due to the fact that the flows were increased by a significant amount from the original model to 
provide a detectable signal, as discussed earlier. Computing the new W (column 3) to conserve 
mass improves this (column 5), but requires strong radial flows, with an rms value about 30% larger 
than the original value (Table [H row 3). Here, the typical region in the convection zone lost its 
mass over a few tens of thousands of years. 


“MF3” is the inverted profile of lHartlep et al.l (j2ni3l ) shown earlier in Fig. [2j The 2nd column 
of Fig. [5] is empty as the radial flow is not inverted for. Thus one expects the mass conservation 
to be poor, and the 4th column verifies this. Interestingly, it’s at the same order-of-magnitude as 
the original profile itself, even withouta Vr- The new Wj shown in the 3rd column, serves to help 
satisfy continuity rather well (column 5). 


“MF4” and “MF5” are the GONG and HMI inversion results, respectively, shown earlier in 
Fig. m The radial flows needed to satisfy continuity at a signihcant level are strong, with maximal 
values in concentrated regions of about 100 ms“^, yet with much smaller rms values. They also 
display a complex pattern of downflows to upflows in latitude, that differ quite a bit from each 
other. In terms of the draining/filling of the convection zone, both inverted solutions are on the 
order of 10^ year. It is somewhat surprising that even with no radial flows the inverted meridional 
profiles from GONG and HMI conserve mass a bit better than the simulated data, but are still 
about 4 orders-of-magnitude worse than the fiducial case. 


What all of this is suggesting is that the flows we are recovering in our inversions can be 
interpreted as being unphysical if our basic understanding of the solar convection zone is correct. 
Stated differently, the retrieved circulation cannot be the whole picture, since we find that the 
radial flows (which were are not inverted for) necessary to satisfy continuity are significant. They 
also differ dramatically from the “standard” picture as in MFl. If these radial flows do exist as 
required to conserve mass, then they should be measureable. Furthermore, the flows predict that 
the convection zone would be drained of mass on much shorter timescales than solar evolution, 
which we are quite certain is not the case. While it is difficult to put too much emphasis on the 
precise numbers here due to numerical inacuracies, in a relative sense this simple exercise suggests 
that the helioseismic results should be carefully explored further. 
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6. Conclusions 


Inversions of GONG S-N travel-time differences that have been ad-hoc corrected for an un¬ 
known center-to-limb systematic effect show the canonical surface poleward meridional flow, below 
which exists a shallow equatorward component beginning at about 65 Mm beneath the surface. 
This is similar to what is seen from HMI data, yet disagreement appears deeper down when a sec¬ 
ond poleward component is found only from HMI. These results may indicate, somewhat curiously, 
that at some level a similar GTL systematic is present in both the ground-based and space-based 
data, which is not an obvious conclusion a priori. 


Irregardless of systematics, convection-zone mass-conservation arguments appear to demon¬ 
strate that the inverted flows are problematic, presumably in both magnitude and structure. This 
is difficult to reconcile, and given the similarities in GONG and HMI results, one obvious question 
is why is there relatively good agreement if they are both incorrect at some level? We note that 
other recent studies using local heliose ismology for nieasur i ng smaller-scale flows have b een more 
challenging than orginally anticipat ed (IDeGrave et al.ll2ni4l : lHanasogell2ni4l : ISvandall2f)15l ). though 
there are promising directions (e.g. IWoodardll2ni4l ). 


Other improvements are possible, such as the development and use of finite-frequency sensi¬ 
tivity kernels, which could also provide a way to potentially measure the radial contribution of the 
circulation. Noise covariances could be more accurately determined, and each day there are more 
and more data that can be analyzed to help reduce uncertainties. Sources of systematics are also 
being explored. Unlike rotation, meridional flows are expected to be a much weaker signal and it is 
not surprising that they have been so difficult to measure accurately despite very hard work in the 
past several decades. Thus, as demonstrated here it is critical to have independent measurements 
from different groups to corroborate any inferred circulation patterns. 


We gratefully thank J. Zhao for providing his inversion results from HMI data, T. Hartlep for 
providing his meridional flow simulation data, and an anonymous referee for an excellent numerical 
suggestion. This material is based upon work supported by the National Science Foundation 
under Grant Number 1351311. A.S. acknowledges support from grant FA-F02-F028 of the Uzbek 
Academy of Sciences. S.K. was supported by NASA’s Heliophysics Grand Challenges Research 
grant 13-GCR1-2-0036. This work utilized data obtained by the Global Oscillation Network Group 
(GONG) program, managed by the National Solar Observatory, which is operated by AURA, 
Inc. under a cooperative agreement with the National Science Foundation. 
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Fig. 5.— Tests of continuity equation for the 5 meridional flow models (MF^) described in the 
text. The columns show the latitudinal speed (column 1), the radial speeds (columns 2-3), and the 
continuity computation (columns 4-5). In each panel the dashed refe rence line encloses th e ±60° 
and (0.7 — 1.0)i?Q interior region. The large flows associated with the lHartleo et al.l (j2013l ) model 
(MF2 and MF3) have been scaled down for plotting to achieve similar values of the other flows (see 
Tabled] for unsealed numbers). 
















- 18 - 


Table 1: Continuity parameters for five meridional flow profiles. 


Row# 

Quantity 

Unit 

MFl 

MF2 

MF3 

MF4 

MF5 

1 

rms 

ms“^ 

1.79e-h01 

1.87e-h02 

1.75e-h02 

9.45e -h 00 

6.77e-h00 

2 

..rms 

Vr 

ms“^ 

3.10e-h00 

1.43e-f 01 

O.OOe-hOO 

O.OOe-hOO 

O.OOe-f 00 

3 

Tyrms 

Ur 

ms“^ 

3.10e-h00 

1.84e-f 01 

2.81e-h01 

4.24e -h 00 

2.99e-f 00 

4 

div pv 0 

- 

3.30e - 05 

2.32e - 04 

4.75e - 04 

-1.73e-05 

2.82e-05 

5 

Y] div pvr 

- 

-3.30e-05 

1.02e- 04 

O.OOe-hOO 

O.OOe-hOO 

O.OOe-f 00 

6 

X](div pve + div pvr) 

- 

-2.63e-09 

3.34e - 04 

4.75e - 04 

-1.73e-05 

2.82e-05 

7 

^ div pVr 

- 

-3.30e-05 

-2.32e-04 

-4.75e-04 

1.73e-05 

-2.82e-05 

8 

yddiv pve + div pVr) 

- 

2.23e-09 

3.69e- 08 

4.90e - 08 

1.83e-09 

2.49e - 09 

9 

y; (div -l-div pvr ) 

- 

-7.97e-05 

1.44e-h00 

l.OOe-hOO 

l.OOe-hOO 

l.OOe-hOO 

10 

(div pve -l-div pvr) 
div pve 

- 

6.77e - 05 

1.59e-04 

1.03e-04 

-1.06e-04 

8.81e-05 

11 

{-p{dp/dt)-^) 

day 

-3.36e-h06 

2.51e-f 01 

3.05e-h01 

-1.82e-h02 

2.89e-f 02 

12 

{-p(5p/at)“Qscaled 

year 

-4.57e-h09 

3.41e-h04 

4.16e-h04 

-2.48e-h05 

3.93e-h05 


Note. — All quantities are computed in the truncated radial region r = 0.75 — 0.995 i?©. The quantities in rows 
11-12 are computed from the divergence of the mass flux using Vr, not Vr, and (...) represents the median of the 
values in the truncated region. Physical units are only provided for relevant quantities. 









